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In each case the time of charging was two minutes. 

The results are given in terms of the total charge in electromagnetic units 
contained in the specimen when the original difference of potential between 
the surfaces was unity. In the first column times measured from the instant 
at which both surfaces are connected to the earth are given. 

These results are plotted as figs. 20-22. For these specimens the values 
of a are, (after two minutes charging) : — 

Glass (test-tube). Paraffin. Ebonite. 

0-0115 0-0048 0-0055 

The author's thanks are due to the Managers of the Boyal Institution for 
placing tbe resources of the Davy Faraday Laboratory at his disposal, and to 
the Director, Prof. Sir James Dewar, for much helpful kindness and 
consideration. 



On Deep-Sea Water Waves caused by a Local Disturbance on or 

beneath the Surface. 

By K. Teeazawa, Ei-Gakushi in the Imperial University of Tokyo, 
(Communicated by Sir Joseph Larmor, F.E.S. Eeceived June 28, 1915.) 

The classical problem of waves produced in deep-sea water by a local dis- 
turbance of the free surface has been investigated by Prof. H. Lamb* in a very 
able manner. He completed the theory of wave propagation in one dimension 
and in two horizontal dimensions when the initial disturbance is concentrated 
in the immediate neighbourhood of a line or a point, assuming that Fourier's 
double integral theorem can be applied in such a case. In this paper I 
venture to discuss the same problem, especial]y the oscillations at the centre 
of the disturbance, following Prof. Lamb's method, in cases where the initial 
disturbance is spread over a certain extent of the free surface ; and, as an 
application of the general solution, I propose to treat the case in which the 
initial disturbing source is situated at a finite depth from the free surface, i.e. 
where the surface wave is produced by an explosion like that of a mine 
under water. 

§ 1. Supposing the water to extend to infinity, horizontally and downwards, 
and taking the axes of x and y on the undisturbed free surface and that of z 

* 'Lond. Math. Soc. Proc., 5 vol. 2 (2), p. 371 (1904), or 'Hydrodynamics,' third edition, 
§§ 236-239, 251-252. 
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vertically downwards, we have, on the assumption that the motion is 
infinitely small and irrotational, 

p J p == d(j>fdt-tgz, (1) 

where p denotes the pressure, p the density of water, g the constant of gravity 
and cj> the velocity-potential which satisfies the equation 

y 2 <£ = 0. (2) 

If % denotes the depression of water-level at the point (x, y, 0) and at time 
t below the undisturbed surface, then the pressure condition to be satisfied at 
the free surface, supposing f to be small, is 



K 



l 

9 



d<f> 



(3) 

z =? o 



and the kinematical condition is 

3f __ r3$" 



dt 
Hence, for z = 0, we must have 



dz Jz = o* 



(4) 



or, if the time factor be e i ^ &t+e \ 



d 2 d> d<b 

dt 2 dz' 

a 2 (b = —g~. (5) 

The solution of equation (2), which makes <j> = when s = oc, is, in 
cylindrical co-ordinates (w, 0, 2), 

^ = «"**J m (ft n r)S}^ ( 6 > 

where J m is BesseFs function of mth order and * is any positive number. 
From equations (5) and (6) we get 

a 2 = #&. (7) 

The modes of irrotational motion of water by surface disturbance present 
themselves, as Prof. Lamb points out, in two forms : (1) by an initial displace- 
ment of the surface, without initial velocity ; (2) by an initial impulse applied 
on the surface, without initial surface displacement. The most general case 
is a combination of these two. 

§ 2. As the typical solution for the case of initial rest we take 

<£ = g !H^ e- k *J m (hsi)GOsme 9 (8) 

£ = — cosa^J w (/m)cosm<9. (9) 
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The water is supposed to be unlimited in extent, and hence there is no 
limitation to the length of waves, but the motion will be the result of the 
superposition of waves of infinite variety of length. Therefore to obtain the 
general expression which embraces the superposition of all such solutions we 
must make use of the double integral theorem of Fourier's type 
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J m (km)MJc f(ot)J m (kct)ada, 
o Jo 



(10) 



with some conditions concerning the function f(m). 
Thus, corresponding to the initial conditions, 

£ = — f(m) GOB 1710, <£=:0, 

we have 
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cj)=:gGOSm0 <&(/<;) e kz J m (km)kdk, 
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f<x> 



■GOB 7)10 <&(k) GOB crt J m (]c&)kdk, 
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where 



'CO 



<&(k) = f(a)J m (hot) adcc ; 



(ii) 



(12) 



(13) 



(14) 



-0 



for these expressions clearly satisfy the conditions specified above. As to 
the disturbing function /(w), we can assume its form in several ways, under 
the condition that the integral theorem (10) must remain valid. 

Now we consider the case where the initial elevation of the surface is 
spread over the whole surface according to the law 

A 



/(«0 = 



m 



(15) 



This is the simplest form which satisfies the double integral theorem. By 
taking the value of A sufficiently small, we may consider that it approaches, 
when m = 0, to the case where the initial displacement of great amount is 
confined to the neighbourhood of the origin.* 

[* It lias been kindly pointed out to me by a referee that this statement requires 
correction and precision. However small A may be taken, the initial elevation is spread 
uniformly along the radius, but is not concentrated in the sense of Cauehy and Poisson, 
though it is to be expected that only the central part when the height is considerable 
can contribute sensibly to the amplitudes of the waves. It might appear that the total 



initial elevated volume, which is expressed by the integral %tt 



fco 



f(&)&d& when m = 0, 



should be finite. But it seems to me that this condition is not a necessary one for the 
validity of the integral theorem (10). This remark applies to both the initial dis- 
tribution A/S5 and A/(6 2 + CJ 2 ) ¥ that is treated in § 4. But for the initial distribution 
A6/(6 2 + S 2 ) 5 that is treated on p. 69 there can be no doubt, for it is of finite total amount ; 

F 2 



Mr. K. Terazawa. On Deep-Sea Water 

Put the expression (15) into (14), then it appears, by a well known 
formula, that 

(16) 



cf> (k) = A J m (ka)da = A/& * 
Jo 



Whence the generalised expressions for <fr and f will become 



(fy = A# cos m# j — -— • e~ ks J m (Jcm)d7c, 

o 



cr 



£ = —Acosm# 
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cos at e kz J m (km) dk 



z — 



(17) 



(18) 



expanding sin at and cos at into series, we have 



00 



(fit 2 ) 



n r» 



<f> = Agt ooa m0 2 ( — 1) w /^xttj- k n e- kz J m (km)dk, 



so 



? = -A<josm0 S(-l)" ( -|^- 

w = ^ • 



w r 



'oo 



k n e~ kz J m (km)dk 
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2 = 



(19) 



(20) 



here <r is eliminated by the relation a 2 = #&. 

The definite integral contained in (19) and (20) can be evaluated as 
follows : — - 

By using Hansen's form of Bessel's function 



J m (kw) 
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e-^wcos^ cos mcjydcj), 



and the integral formula 



co 



k n 6 ~ k ( z + ^ cos $) dk = 



?i! 



o 



(# + iw cos <£) 



»+] 



, (»>0) 



we find that 



If we put 



then we have 



k n e~ kz J m (km)dk = r— r .^ B+r 

V y 7T- J (3-f»COS<£) n + 1 



r = v / (^ 2 + ^ 2 ) J ^6 ~ s/r, 



(21) 
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k n e~ kz J m (kxs)dk = 



^ ! 



7T? 



.w+1 



7T 



cos mc£ <:/<£ 



fri, + 0i, 3 -l)*cos£l n+1. 



and when 5 vanishes it is concentrated at a point, the result being then the well known 
form of Cauehy and Poisson. 

A main object of this paper is, as stated, to trace the effect on the waves propagated 
to a distance, and especially on the course of the disturbance at the centre, of a 
broadening of the initial disturbance from a mere point to a wider region. This has 
also the advantage that we do not have to proceed to a limit.] 

* Gray and Mathews, 4 Treatise on Bessel Functions,' p. 79. 
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Now, define the function V n m {fJb) by 



n 



p„-»0*) = (i -^r m/2 . ..v % {p)dfi m , 



fX 



P n (fi) being the zonal harmonic of nth order, then the integral expression of 



P —7)1 

L n 



(fi) } which corresponds to Laplace's integral for V n (/a), will be 



P»-"V) 



i m n I 



7T 



cos m<f> d<f> 



7r(n + m)\} [/a + (^ 2 — -1)*cos<£] w+1 



Therefore the required integral becomes 



■°k-**j m (& w )dfc = ( -^|^ p„-»o*) 



(22) 



The result of this integral, in the case of n = m, was found by Prof. E. W. 
Hobson.* In this case, if we remember that 

P.--0*) = p=^, P,»(<*), 

(n + m) ! 
where P w w (^) is the associated function defined by 

(XLL 



we have 



??z) 



(23) 



in agreement with his result. Hereafter we will use (22) for n^m and (23) 
for n~m. 

Hence, putting these values (22) and (23) into equations (19) and (20), 
we get 

>r ( } (2*M-l)lW/ * (M) 



* 



r 



(2w+l) 



and 



^-^^C^-D-'-^^r^) 



^l-'^v^ 



(24) 



(25) 



?i = HI 



* * Lond. Math. Soc. Proc.,' vol. 25, p. 73. 
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The values of V n m (0) and P w m (0) are ( as follows : — 

(m + w) " 7 



■^ 



mi 



i' 



p m /Q\ __. / ^ \ (»- m)/2 j; ; ^M.(^ + m-l) 

2 .4...(w — m) ; 
= 0, 



(w=0) 

(%—meven) 
(n—m odd) 



J- (26) 



= 1 . 3...(2?i— -1), 

For the case of symmetry about the origin, i.e. m = 0, the first summations 
in (24) and (25) disappear, and therefore we have simply 
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m — 



r n to K } (2rc+l)!\r 



1\(V>), 



fcm = 



w L 4 ! \ w / 8 ! \ to- 



(27) 
(28) 



similar expressions to those which were found by Prof. Lamb for the case 
of complete concentration of the source. 

The next simplest case is that of m = 1, though it may not be very 
important in practice. In this case we obtain 



* 



m = 1 



Agt cos 8 



p ^> + 2 (-^^(t^ 



n=l 



n(p)j> 



(29) 



b m = 1 — " 



A cos 6 



m 



2!W/ 6! \w/ 10! \W " J' 



(30) 

These solutions do not give any information as to what takes place at the 
immediate neighbourhood of the origin. Not only the initial data but the 
displacement at any time will be infinite at the origin; such a point is 
excluded in our fundamental assumption. The appropriate initial data and 
the solutions to illustrate this point will be given presently (§ 4). 

§ 3. The series (25) is not convenient when we deal with the case in 
which gt 2 /m has large values, since it converges rather slowly in such a case. 
The suitable expressions in this case, at least when m — and m = 1, can be 
obtained by a similar method to that employed by Prof. Lamb for the case of 
complete concentration of the source. 

If we put 

T m (z)'= cos aie~ kz 3 m (km)clh, (31) 



then returning to (18) we have 

f = — Aeosm0 . T m (0). 



32) 
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By using the integral expression for J m 

1 [" 

J m (Jem) = - cos (m<l) — Jem sin <f>)d<f>, 



and putting 



•co 



G(z) = I e~ kz cos at cos (Jem sin <p)dk, 

o 



■oo 



B(z) — e~ kz cos at sin (Jem sin $)c&fc, 



o 



with 
we get 



a 2 = gk ; 



T m (z) = —I C(V)cosm<jf><^>4- - S(»)sinm^d^. 



Now expanding cos at into a series and using the formulae 
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Jc n e~ kz cos xJcdJe 



n\ 
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P 



n+ 



jCOs(n+l)yfr 
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h n e~ kz sin xlcdk = — n sm(n + Y)y}r 



o 



_p 



> 



where 

and 
we find 



p =z x /(z 2 -{-x 2 ) i ^r — tan l (x/z) 
x = OTsin^>; 



3 



00 t? ! int 2 \ n 

s(2)= 3o ( - i)B 2 ;^(7T sin(M+i) ^ 

= corresponds to p — m sin <£, t|t = Jtt. Hence 



O(0) 



CO 1 
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sin0 
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x.o.O jl , o . d . i . y 



S(0)=-4--| 
where co = gt 2 /2msmcj). 
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(33) 



(34) 



(35) 



The functions C(0) and S(0) are closely connected with FresneFs integral 
as is well known. If we put 

1 „, N 1 



C(0) = 



S(0) = 



m sin cjy 
m sin<£ 



M(©) 



-j 1 — N(g>) 



^ 



(36) 



> 
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then we shall have' 



M(o,) = ^{cos| 



m u du . . &) f w . u die 
cos - —-T- + sm - sm - 



■^ 



N" (w) = ^~i sin ^i eos 
2 L 2J o 



o * v ^ 

w % du 
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^ J 



2 <\/% 
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COS- 
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2 y/% 



^ (37) 
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If gt 2 j2i$ is a large quantity, then &> is also large, and the functions M 
and N" can be expressed as semi-convergent series of the formf 



M(ft>) = v v 2 ; - |sin - +cos - 



1 1.3.5 



ft> 
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If we neglect a> * and higher powers of it in these series, we have, from 
36), approximately 



O(0) 



2m sin (ft 
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Putting 
and 
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_sin <j>) J (sin $) 3/2 
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T m (0) can be written in the form 



T«(0)=i^*{U in + V«}. 



(41) 



From (40) it appears that U m = when m is an odd number, 
and V m = when m is an even number. 

If we find the values of U TO and Y m in proper form, we can obtain the 
expression for f by (41) and (32). 

The special case m = was discussed by Prof. Lamb and it was found 
approximately that 



T„(0) = Aos(^\ 

US \4rT557 



(42) 



* JLamb, i Hydrodynamics/ p. 366. 

t Lord Eayleigh, ' Scientific Papers,' III, p. 129. 
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and therefore we get 



£ m = 0= -V^Aoos^Y (43) 

For the ease m = 1, TJi = 0, 



Vi = 



,?l * ' ' OL \ . I OL 



d(j> 



cos i -r— r — sm . • ,, r , . , , 
sm 9 / \sin<£/J -y/SlIKp 





\ -1 



and TxCO) = -f^X (44) 

The value of the integral contained in the expression of Vi cannot be found 
very easily. Since we see that 

,7r sin / g V # = 2 ftin f » \ _*£_ = 2 T ^^ ^ 
. \sin </>/ v/sin <£ J \cos </>/ ycos <£ J ± </u ^/ (u 2 — 1 )* 

and similarly 



ff ' a \ deb aT cos ait clu 



o 



cos ( - — - j , / , = 2 i 



v sin <j)J ^/sin <\> Jj y'-w, y^ 2 — 1) 



and if we remember that 

sin x /tt T / \ cos x /ir T / \ 

^f = V 2 • J * (a!) ' ^jr = V 2- J - (a;) ' 

then we can evaluate these integrals, following the method which is fully 
explained by 1ST. Melsen in his book on Cylinder-Functions.* I calculated 
these integrals after this method, and found that 



K45) 



r^ v^-l) =-«*W2) ,V8 {j i 8 («/2)-v / 2 J i («/2)J_ i (a/2)} ' 

" ^T^TI) = «*(W2) 3/2 {j_, 2 («/2)- v /2J i («/2)J_ i («/2) j 

Eor large values of a we can use the asymptotic expansions of Bessel's 
functions, and, to the degree of approximation above mentioned, it is 
sufficient to retain only the first term of them. Thus we may take 

Jr (a/2) = ,J - cos (*/2-3tt/8), J .(a/2) = --I- cos(a/2-7r/8). 

Putting these values in the above formulae (45), then by (44), after 
simplifying, we have 

Ti(0) = - v/2/w . sin(^ 2 /4w). (46) 

Substituting this for Ti(0) in (32), we obtain 

f m = 1 = v /2A/w . sin (gt 2 /4:m) cos 8. (47) 

* ' Handbuch d. Cylinderfunktionen,' §§ 77, 81, 82. 
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Though we have been considering in the above only the terms involving 
cos m#, it must be understood that the same analysis could be applied to the 
terms containing sin md, if the initial data be so given. 

§ 4. As a proper assumption to illustrate what occurs at the origin, take 
the disturbing function of such a form as 

/( ro ) = A/v / (^ + ^ 2 ) (48) 

instead of (15). This function does not violate the double integral theorem 
(10). Supposing the value of h to be sufficiently small, this approaches, at a 
distance from the origin, to the function assigned in § 2, and yet it is finite 
at the origin. Therefore we might take, even in this modified assumption, 
the solutions in § 2 as those for the point at which w has a moderately 
large value. 

For the case of symmetry (m = 0),* the function <I> (k) defined by (14) 
becomes 

This integral can be effected by using the double integral theorem (10) in 
a special case m = 0, and is found easily to be 

®(k) = Ae~ bk /lc,f (49) 

supposing h to be positive. 

Hence we have the expressions for $ and £ from (17), (18) 

4> = Ag P* S ~^ e-<* + WJo(lsm)dlc, (50) 

o a 



"00 



f = — A cos at e~ bk J (lcm)dh (51) 

Jo 

Making use of the integral formula (23), we get the following solutions, in 
place of (27) and (28) 

^^(-^(^(f)"^ (52) 

with s = x /[(z + If + m 2 l v = (z + b)/s, 

and £ = - - £ (-l) n ^(^- a ) Wp »(^ ^ 53 > 

with s = -v/C^ + ro 2 ), vq = 6/s . 

At the origin (w = 0), we have 

S = ^ z; = 1, 

* The integral theorem (10) is valid f or So — only when m = 0. 

t Or, as a special case of the integral of N. Sonine ('Math. Ann.,' vol. 16, p. 50), we 
get the above result, putting m = 0, n = - J in his formula « 14 . 
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and since V n (l) = 1, the equation (53) gives 



2\n 



If we put, for the present, 

gt 2 /2b = x, 

then the above expression can, referring to (35), (36), be put in the form 

£o = A/5 . {-l-iM^ + N^)}, 

and, therefore, using the integral expressions (37) for the functions M and 1ST, 
we have 

£ = A/h . | — 1 + |\A e ~* /2 eFl*dx/</x j . 

Or, putting again, for the sake of simplicity, 

Ja? = t 2 = ^ 2 /4&, Jb = Ay/b, (55) 

it follows that 



7] = — \-\-2re 



~ T 2 



e r \lr. (56) 



o 



The expressions (54) and (56) for the displacement of the point which 
was initially situated at the origin are equally inconvenient to discuss its 
subsequent motion, but it is not very difficult to obtain the general feature 
of it from the latter form. For the small value of r the value of the 

"T 

function e~ 72 e T2 dr can be obtained by expanding the exponential 

Jo 

function in the integral sign and integrating it term by term. For the 

large value of r the asymptotic expression for the function e~- r *\ e r2 dr is 

Jo 

1 1 1 3 '1 3 5 

±_4- _ -{- L ' J- 1,0,,J l- * (^X 

2t 2 3 t 8+ 2V + 2V + "" K } 

Therefore the course of the ordinate 7) at the origin tends approximately to 
the form 

1 46 , 1 . 3 /46\ a , 1.3.5 ^ ox 

* = 2^ + *-w + -^-y + - (t>8) 

when £ is great. 



* E. T. Whittaker, 'Modern Analysis, 5 p. 170; A. E. H. Love, 'Phil. Trans.,' A, 
vol. 207 (1907), p. 195. 
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Table I. 





f T 




*T 




/*-M0 


T. 


e—r 2 


&*dr. 


l — 2re~r 2 


er 2 dr. 


i- + (l-2r 2 )e-T 2 


er % dr. 




w 



00000 


Jo 
+ 1 -ooooo 


*> 





o-o 


o- 


o -ooooo 




o-i 


-09933 


-98013 


+ -19744 


0'2 


-19475 


-92210 


-37917 


0-3 


-28263 


-83042 


-53176 


0-4 


-35995 


-71204 


-64477 


0*5 


-42443 


-57557 


-71222 


0-6 


'47476 


-43029 


-73293 


7 


-51034 


-28552 


-71021 


0'8 


-53210 


-14864 


-65101 l 


0*9 


-54073 


+ 0-02669 


-56475 j 


1-0 


-53808 


-0-07616 


-46192 


a/2 


-4525 


-2799 


+ 0*0567 


a/3 


-3640 


-2609 


-0-0879 


2 


-3013 


-2052 


-1091 


a/5 


-2585 


-1560 


-0904 


a/7 


-2075 


-0980 


-0517 


3 


-1783 


-0698 


0-031 


5 


-1021* 


-0213f 


-0045 X 


10 


-0503* 


-0 '0051f 


-0*0005 % 



# Calculated by using (57). f Calculated by using (58). J Calculated by using (86). 

By this Table the curve representing approximately the displacement at 
the origin can be easily traced, the unit of time and displacement being 
modified by (55). 




Waves caused by a Local Disturbance. 69 

The highest point in the curve corresponds to the initial prescribed eleva- 
tion at the origin. As r increases from zero the elevation decreases and 
reaches once to the value nil ; and then begins the depression of the surface 
at the point under consideration. After that there is only one maximum of 
the depression, the amount of which is smaller than that of the initial eleva- 
tion. Then it decreases more and more, very slowly, until, after an infinitely 
long time, it takes the limiting value zero. The number of zero-point and 
maximum-point can be determined without using the above figure. 

We can assign several assumptions to the disturbing function to illustrate 
the movement at the origin. For example, if we take the form 

f(m)^Ab/(h*+m 2 )W (59) 

instead of (48), then it appears that 

<£(&) = Ae~ bk , (60) 

which can be found in a similar way as before. 
In this case 



f = - A cos o-te~ M J (Jm)kdJc, 
Jo 



and, performing the integration, we have 



^ = - l! 1 ( - irl (-2^Wi©"" lp -^* < 6i > 



where s and v are defined in (53). 
At the origin this becomes 

* A 3 / -,n„ i n\ /gt 2 \ n - 1 

Concerning this series and also as to the motion at the origin, we can make 
a similar interpretation to the former case. 

II. 

. § 5. Quite analogous treatments can be employed in the case where the 
initial condition is an impulse applied on the free surface. The typical 
solution is 

pcf) = cos at e~ kz J m (lm) cos m0, (63) 

£= ljgp . <r&m<rtj m (km)oo&m9. (64) 

[* This formula gives, in fact, the Cauchy-Poisson series as a first approximation 
when b is taken very small (Lamb's ' Hydrodynamics,' p. 408, equation 20). See the 
footnote on p. 59.] 



70 Mr. K. Terazawa, On Deep-Sea Water 

Corresponding to the initial conditions 

pcf) = F (m) cos m6, % = (65) 

at the free surface, the general solution is 

cos md rco 



<b 



P 



V(k) cos <rte~ k *J m (km)k elk, (66) 

o 



•oo 



^ = cosm^ y^ a gin ^ 5m Q. m ^ mh . ( 67 ) 

#P Jo 

*QO 

where M^fc) = F(a)J TO (Aja)otrfa. (68) 

Jo 

At first we assume the initial impulse is given by 

F(w) = A/w, (69) 

then the value of M/ 1 (/c) becomes 

^F (ft) = A/&, 

corresponding to (16) in the former case. 
The explicit form for the solution can be obtained by performing the 

operation -—-57 upon that of the former case. Thus from (24) and (25) 
99 ot 

we get 

, Acosm# f^ 1 / i\„(^ + w)!/^ 2 \» •, x 



+ |(-1)» (^^^"p^^)!, (70) 



n—m 



where r = \/(z 2 -f- ^ 2 ), ^ = s/r ; 

£ = Acosmg f ^z jxn-x (^+m)! /^ 2 \ n p w - ffl , v 
p^sr L»=i (2%~-l)lW/ 



+ 1 (- 1 ) n " 1 (^E§t($)" p » ,,l ( ) ^ (71) 



n=m 



Specially for m = 0, from (27) and (28), 

^-ph-^ffi™' (72) 

y A fl a (gP\ a 1 3 .3 3 /W, 1* . S* . & (ffiW 1 m . 



and for m == 1, from (29) and (30), 

A ^{po-v) + i<-i)*M 1 t- x " 



^ = *<^ Po - V)+ s(-l)»i^i ^ P.W K (74) 



G^ = A^i {l».^»' + ^p , ....l (75) 
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Now, when gt 2 jm is a large quantity, we can make use of the function T 
defined by (31). Thus 

f — — ^ cosm 0dT m (O) 
pg dt 

Eef erring to the expressions for T (0) and Ti(0) found in (42) and (46), we 

get 

Km = o = /0 y » sin (gt 2 /4&) (76) 



for m == ; and 

\/2pw 
for m = 1. 



&» = 1 = /o__2 C0S C^V 4t;T ) C0S # C 7 ?) 



| 6. As in § 4, we take the function expressing the initial impulse of the 
form 

F(*r) = A/v/(& 2 + ^ 2 ) (78) 

to illustrate the history of the point which was situated initially at the 
origin. 

For the case of symmetry about the origin (m = 0), the function M^'/fc) 
becomes 

^F (k) = Ae-**/ft. 

Thus we have 

Jo 



' /> 



CO 

cos <r* g-(* +6 >* Jo(Aot)^, (79) 



f = — f W ff sin <r* . e™ 6 * Jo (few) & (80) 

Or, performing the operation — 37 upon (52) and (53), we obtain 

gp at 

*=%t!„<- 1) "ln(7) p -M- (si) 



in which 



s = ^[(3 + b) 2 + m 2 l v = (z + 6)/s, 
% = \/(b 2 +m 2 ), vq = &/s 



as before. 

At the origin (w = 0), f becomes 



go = ^7 S(-l)»-* • , ^ . (83) 
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This series can be put in a similar form to that in § 4 ; or rather directly 
from (56) it follows that 



fT 



v = t + (1-2t 2 >- t2 efdr, (84) 



o 



in which, for simplicity, it is put 

v = (pb/A)^/(gb)U r = t^ifffib). (85) 

For large values of t, using (57), we can take for 77 its asymptotic expansion 
of the form 

1 2.1.3 3.1.3.5 4.1.3.5.7 



2t 3 2V ' 2V " 2 4 t 9 



(86) 



From these expressions we can make out the general feature of the move- 
ment of the point under consideration. 

By the aid of the. Table I we can draw the approximate diagram which 
represents the displacement at the origin, the units of scales being modified 
by (85). 




For t = 0, we have rj = and drj/dr = 4. Therefore as soon as the impulse 
is applied on the free surface the point at the origin gains suddenly a finite 
velocity and begins to move downwards. It will be seen, from the above 
figure, that at about r = 0*6 the depression becomes maximum, at about 
t-— 1*5 it takes the value zero and then begins the elevation of the surface, 
and it takes the maximum value at about t = 2, after that it decreases 
slowly until it tends asymptotically to its limiting value zero, and the point 

comes to rest. 

If we assign another assumption to the initial impulse like that expressed 
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in (59), it will probably take place by a similar mode of motion, but in a 
somewhat complicated form. 

III. 

§ 7. Now we come to the problem of surface waves caused by an explosion 
which takes place at a finite depth from the free surface.* In this case, we 
might take the value of the impulse, i.e. that of p(j>, given by the explosion as 
the initial datum, but the mathematics would be complicated and it might 
not give any concrete results. The following consideration, however, may 
serve to treat the problem without touching on the difficulties. 

Suppose that an explosion under water sends out a pressural pulse of the 
type of a wave of sound, propagated by compression at high speed, and 
leaving unruffled water in its wake. As each ray reaches the free surface it 
is reflected totally, as required by the condition that the impulse is null at 
the surface, while a surface normal velocity caused by the pulse remains and 
determines a hydrodynamic flow throughout the water.f We may take the 
velocity of the propagation of the pressural pulse to be infinite, then the 
image method in electrostatics will be applicable as follows. 

Let the point at which the explosion occurs be at the depth h from the 

free surface and on the axis of z. Suppose the initial value of <£ due to the 

source at C to be A/?*, where r is the distance of any point in the water from 

C, then we must consider a fictitious negative source of equal strength at the 

image of C with regard to the free surface, to fulfil the condition that the 

initial impulse is zero at the free surface. The initial value of $ should 

then be 

= A(l/r-l//), 

where r = \/[(z - h) 2 -f w 2 ], r f = \/[(z -f h) 2 + m 2 ~\. 

This value of cj> gives the initial normal velocity 

-3<£/&= -2Ah/(h 2 + m 2 yi 2 (87) 

over the free surface and it dies away of itself. Thus the problem is 
transformed to one of determining the surface motion, being given the initial 
surface normal velocity (87). 

If we take 

<£ = cos at e~~ kz Jo(Jm) 

* This problem and the next following were undertaken mainly to illustrate in a 
general way the nature of the tidal waves produced by a submarine earthquake, with its 
source at a point or on a long straight fault in the strata. 

t The source of the explosion is supposed to be so deep or its force so gentle that the 
surface of the water is not broken by the ejection of a column of water. 

VOL. XCII. — A, a 
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as the typical solution, then the general solution corresponding to the initial 
data (87) will be 

f°° f°°r 2A7^ "1 

cfr = cos crte~~ kz Jo(k® )clk\ —th> 57^ Jo(ka)ada (88) 

Jo J L (/r+<n 2 J v > 

with ar = \/(gk). 
Since, by (60), 

.CO 7 

(fta + w 3wg J o(fe t )«rfflt = e~ hk , 
we have $ = — 2 A cos Gte"~^ +h ^J Q (km)dk, (89) 

Jo 





.oo 



2A 
and accordingly £= -I o-sino-^e kh J (km)dk. (90) 

# Jo 

These are the same expressions as (79) and (80) up to the constant factor. 

Therefore putting 

s = ^/[(z + hy + m 2 ], v = (z + h)/s, 

% = a/(^ 2 + ^ 2 ) ^o = A/s , 

we get <f> = " ^ 2 o (-l)»|^(^)*P»(^ (91) 

and' f = -M |(_i r i^^^\"p n ( Fo) . (92) 

At the point just above the explosion, that is at the origin, w = 0, putting 

g ty4:h=T a , (93 



we have 5> = ~^?){ T + ^ ~ 2t2)6 ' T ' \f dr } ' (94) 

Comparing (94) with. (84), we see that the curve representing the displace- 
ment at the origin is quite the same as in fig. (2), only inverted in form and 
the unit of ordinate being different. 

If we take 

Tl — 0*6, r 2 = 1'5, t 3 = 2, 

as the values of t which give the maximum, zero and the minimum of the 
displacement, and g — 32, then the times at which they occur will be 
given by 

respectively, in which tfi, t 2 , and £ 3 are expressed in seconds and h in feet. 
Conversely, if we can measure any of them, then the depth h can be 
calculated approximately by using the above relations. 
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§ 8. We will conclude this paper by discussing the explosion problem in a 
horizontal canal of infinite depth and infinite length. Take the axis of x 
along the canal and on the undisturbed free surface of water, and that of z 
vertical downwards as before. If we suppose that all the circumstances are 
independent of y, we must consider the initial value of <£ given by the 
explosion to be of the form — Alogr.* Therefore, considering the image as 
before, we have the effective initial normal velocity at the free surface 

-c<j>/dz = -2Ah/(h 2 +x 2 ) (95) 

in place of (87). Corresponding to this initial velocity, the general solution 
will be 

(h = — I COS 0^0 kz QQ&KX-=- 





in which a = \/{gk). 

By a well known formula 



o 



h 2 + a 2 „ 



cos Jcada, 



00 h cos lea 7 , ,. . 

-72- 2 tta = *7r6 M ,f 



'00 



we get <£ = — • 2 A cos er£ g~(*+ ^)* cos fe dk/k, 

Jo 

O A Too _ ayr\ rvf 

and consequently f = — = — e~ hk cos kx dk (96) 

^ Jo * c * 

Expanding sin cr£ and then making use of the integral formulae (34), it may 
be shown that 

?=s -^.^-^^^(f)" 008 ^ 4 - 1 ^ (97) 

where q = v/(^ 2 + A 2 ), i/r = ta 1 ii~ 1 (x/h). 

At the point just above the explosion, i.e. along the axis of y, x = 0, 

?0_ A -io ( } (2» + l)lUr 
This series can be put in the form 

5> = - -^ «-** fW-r, (98) 



in which r is defined by (93). 



o 



* Since this work was finished my attention has been called to a paper by Prof. Lamb 
in the l Annali di Matematica,' vol. 21, p. 237 (1913), where the surface waves that 
accompany a cylinder travelling uniformly transverse to its length under water are 
investigated. 

t H. "Weber, ' Die part. Diff.-Gleich,,' vol. 1, p. 43. 

Gr A 
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For large values of t, if we use the asymptotic expansion (57), then 

^-!H i+i (f) +i - 3 (iy +i - 3 - 6 @) ,+ - 

Making use of Table I we can trace the curve for £o. 




In this figure the units of length and time are modified by 

fi> = 4Ai7o/vW> ^ a /4* = **- 

From this figure it will be seen that the line which was initially situated 
along the axis of y is suddenly put in motion and there begins an elevation of 
the free surface, and that after a certain time (about t = 0*9) there is the one 
only maximum of the elevation, and then it decreases more and more until 
the surface takes the original plane form. 

If we take r = 0*9 for the value of r which gives the maximum of the 
elevation, then the time at which it occurs is given by 

t = 045^(1^), 

in which t is reckoned from the time of explosion in seconds and h is measured 

by feet. 

Lastly, at a point very far from the source of disturbance, if we neglect h 
compared with x, we may put q = x and ^ = tt/2 in (97). Then we have 

2Atf 1 (gfi\ 1 /gg\» - 1 _/ggY„ 



^« "" x 11. sW 1.3.5.7 W + 1.3. 5.7.9. ll\2a? 

(99) 

* This formula and those that follow are applicable only when gt 2 /2a; is moderately 
small, this being the only range of values which matters in the physical problem, each 
phase of the physical disturbance travelling, as is well known, with constant acceleration. 
When gt 2 /Z,v is large, the most essential terms of the series (99) are situated far along it, 
and in deducing them from (97) we must not replace cos(w + l)^ by unity. Thus the 
impossible physical feature of persistence of the amplitude for all time with period 
diminishing rapidly without limit, which is indicated in (100) and in fig. 4, would not 
occur in an exact solution : it can be cut away while the earlier undulations of the curve 
remain applicable. 
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17*7 



or referring to (35) and (36) this can be put in the form 

£, = -2At/x . N(©)/a> 

with co = gt 2 j2x. 

For large values of gt 2 /2x } from (38), we get 



£. = -W 



'2tt 



f . off gt 2 \ 4A f' - J2x\* 

< sm^ cos*— v < 1 — 1.3 --5) +... 

L 4^ 4x J gt L \gt 2 / 



* 



(100) 



Further, if we put 

then 

Again, if we put 



ZA 



\/W2)C. <' = Vfc/2*)> 



?? 



00 



N<r 2 )/^ 



(101) 



(102) 



then we have ^ = — N(l/a?'). 

The function JS(z)/z has been tabulated by E. Lommel in his memoir on 
the diffraction problem,! and the other functions needed are deduced in 
Table II, so that the forms of the few waves can be traced without difficulty. 

Table II. 



z. 


3. 


N(«) 

Vz. 


BT(*). 



1 
2 
3 

4 
5 


-oooooo 

+ -323905 
-593492 
0-765194 
-814623 
-741108 


-ooooo 

+ -32391. 
0-83933 
1 -32535 
1 -62925 
1 -65717 


-ooooo 

+ 0*32391 

1 '18698 

2 -29558 

3 -25849 
3 70554 


6 
7 
8 
9 
10 


0-567111 

-333027 

+ 0-088400 

-0-118539 

-250619 


1 -38913 

-88110 

+ -25003 

-0-35562 

-79253 


3 -40267 

2-33L19 

+t) -70720 

-1-06685 

2 -50619 


11 
12 
13 
14 
15 


'288825 

-235250 

-0 '111413 

+ 0-047520 

-201198 


-95793 

-81493 

-0-40171 

+ -17780 

-77924 


3 -17708 

2 -82300 
-1-44837 
+ 0*66528 

3 -01797 



* The first term corresponds to the known Cauchy-Poisson expression for the disturb- 
ance at a distance due to a point source, and the second is the modification in it that is 
produced by the finite extent of the source : it represents a gradual fall of the surface 
to the general level. 

t < Abh. d. k. Bayer. Akad. d. Wiss., 2te CI.,' vol. 15 (1886), Tafel V. 
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Table II — continued. 



z. 



NO) 



z. 



N(«) 


Vz. 


1 -25312 


1 '48142 


1 -40641 


1 -04478 


+ 0-48353 


-0-14134 


-67817 


-99676 


1 -02025 


-74397 


-0-23655 


+ 0*37679 


-94501 


1 -32844 


1 -43138 


1 -22912 


-76983 


+ -16529 


-0 -43709 


-83103 


1 -08444 


-97210 


-58150 


-0-00881 


+ 0*60525 


1-10984 


1 -38093 


1 -35170 


1 '02888 


+ -49105 


-0-13050 


-68400 


1 -03431 


1 -09606 


-85448 


-0-36906 


+ 0-24100 


-82602 


1 -24242 


1 -38795 


1 -22668 


0-79779 


+ 0-20605 


-0*40406 


-88333 



N(*)- 



16 

17 
18 
19 
20 

21 
22 
23 
24 
25 

26 
27 
28 
29 
30 

31 
32 
33 
34 
35 



51 
52 
53 
54 
55 

56 

57 
58 
59 
60 



0-313279 
-359296 
-331493 
-239689 
+ 0-108121 

-0-030844 
-144587 
-207839 

-208258 
-148793 

-0-046392 

+ -072514 

0-178590 

-246628 

-261363 

-220766 

-136088 

+ 0-028774 

-0 -074961 

-150517 



36 


-180741 


37 


-159812 


38 


-094332 


39 


-0-001410 


40 


+ -095698 


41 


-173327 


42 


-213082 


43 


-206133 


44 


-155109 


45 


+ '073202 


46 


-0*019241 


47 


-099771 


48 


-149290 


49 


'156580 


50 


-120841 



-0-051678 

+ 0-033421 

-113462 

-169071 

-187150 

0-163922 

-105670 

+ -027050 

-0-052604 

-114037 



5 -01246 

6 -10803 
5 -96687 
4 -55409 

+ 2 -16242 

-0 -64772 

3 -18091 

4 -78030 

4 -99819 
3 '71982 

-1-24620 
+ 1 -95788 

5 -00052 
1 -14221 

7 -84089 

6 -84375 

4 -35482 
+ 0-94954 
-2-54867 

5 -26810 

6 -50668 

5 -91304 

3 -58462 
-0-05499 
+ 3 -82792 

7 -10641 

8 -94944 

8 -86372 

6 -82480 
+ 3 -29409 

-0-88509 

4 -68924 

7 -16592 
7 -67242 
6 -04205 

-2-63558 

+ 1 -73789 

6 -01349 

9 -12983 
10 -29325 

9 -17963 

6 -02319 

+ 1 -56890 

-3-10364 

6 '84222 



The annexed curves are the graphical representations of equations (101) 
and (102), which are traced by using the above Table. They enable us to find 
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out how the displacement of the surface at a great distance from the source 
of the explosion varies for different times and different places respectively. 
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The figure 4, which belongs to equation (101), shows the rise and fall of the 
■surface at various times and at a particular place of a great distance from the 
origin, the unit of time being modified by tf = t^/(g/2x). For moderately 
rsmall values of t', which alone are applicable, the curve is unsymmetrical with 
regard to the axis of time and as tip time increases the rise and fall becomes 
more and more symmetrical and, on the other hand, the oscillations follow 
one another with ever decreasing period. For larger values of t', the course 
of the curve ceases to represent the motion, which really dies away, for 
reasons explained already (foot note supra). 

As a practical application, the earlier part of this curve might serve to 
determine the distance of the point under observation from the source of the 
explosion, if we measure the time of any zero-point of the displacement or the 
interval (i.e. period) between any two adjacent ones. For example, if we take 
4f = 2*9, 3*7, as the values of t' of the zero-points^ and if we measure the 
actual time 6 between them, then the distance x will be given by 

x = #<9 2 /2(3-7-2-9) 2 = 25<9 2 

in the foot-second-unit system. It will be noticed that the period 6 in this 
formula would be a large quantity of the order ^x/5. 

The figure 5, which belongs to the formula (102), shows the wave profile at 
a particular time at various distances from the origin, the unit of horizontal 
scale being modified by x' — (2x/gi 2 ) x 100. The second part of the diagram is 
condensed along the axis of x f to 1/10 scale of the first part. As we advance 
farther and farther towards the infinite distant point from the origin, the 
essential groups of the waves are found to increase continually in length and 
to diminish continually in height, in such a manner that t 2 jx remains constant 
in each phase. This figure applies only at great distances from the origin, 
i.e. when x is large, so that t has to be large for the earlier part to be applicable ; 
e.g., x' — 2 gives t = l'8y/x in the foot-second-unit system. 

In the region some distance away where equation (102) holds, and therefore 
the fig. 5 is applicable, the wave-length can be determined by finding the 
values of x' of the earlier zero-points in this diagram. 

In conclusion I wish to express my sincere thanks to Prof. Sir J. Larmor 
for his kind advice during the progress of this work. 



